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ABSTRACT 

Previous stressed oxidation tests of C/SiC composites at elevated temperatures (350°C to 
1500°C) and sustained stresses (69 MPa and 172 MPa) have led to the development of a finite 
difference cracked matrix model. The times to failure in the samples suggest oxidation occurred in 
two kinetic regimes defined by the rate controlling mechanisms (i.e. diffusion controlled and 
reaction controlled kinetics). Microstructural analysis revealed preferential oxidation along as- 
fabricated, matrix microcracks and also suggested two regimes of oxidation kinetics dependent on 
the oxidation temperature. Based on experimental results, observation, and theory, a finite 
difference model was developed. The model simulates the diffusion of oxygen into a matrix crack 
bridged by carbon fibers. The model facilitates the study of the relative importance of 
temperature, the reaction rate constant, and the diffusion coefficient on the overall oxidation 
kinetics. 

INTRODUCTION 

Ceramic matrix composites, specifically continuous carbon fibers in a silicon carbide 
matrix, are proposed for high-temperature structural applications due to their high-strength to 
density ratios and ability to withstand high temperatures. High temperature capability allows 
engines to operate more efficiently, and with a reduction in the complex cooling systems often 
required for metallic alloys. Performance also is increased by weight reductions from lighter 
materials. Fiber reinforcement of the ceramic matrix allows the inherently brittle material to fail in 
a more graceful manner due to crack bridging and fiber pull-out. These properties make ceramic 
matrix composites ideal materials for inert environments, however, carbon fiber reinforced ceramic 
matrix composites are prone to corrosion in oxidizing environments. 

Many studies of the oxidation of carbon fibers and pyrocarbon interphases in ceramic 
matrix composites involve thermogravimetric analysis (TGA) in which weight loss of the 
composite is monitored over time [1, 2, 3, 4], These types of studies can provide valuable 
information. However, the oxidation of the carbon constituents hi the composite for unstressed 
and stressed states is very different. In an unstressed state, silica forms near 1200°C and cracks 
close near the processing temperature. These effects protect the interior of the composite from the 
outside environment [2, 4, 5], However, in real application conditions, wide temperature ranges 
and stresses will prevent cracks from sealing. 

The oxidation of carbon fibers in a SiC matrix has been studied experimentally hi the 
stressed state [6], Temperatures ranged from 350°C to 1500°C and stresses of 69 MPa (10 ksi) or 
172 MPa (25 ksi) were applied. Although weight could not be monitored as in TGA experiments, 
the times and strains to failure of the composites provided an indication of the extent of oxidation. 
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Times to failure suggested two regions of oxidation kinetics with a strong dependence on 
temperature. In the temperature range of 350°C-750°C, times to failure dropped sharply as 
temperature was increased. Samples went from surviving beyond 25 hr. with a high residual 
strength to failing in 91 minutes at 69 MPa and hi 22 minutes at the higher stress of 172 MPa for 
temperatures of 750°C. In the higher temperature range of 1000°C-1500°C, temperature 
differences did not have as significant an effect on sample life. Sample lives remained near the 
same level as those at 750°C with some deviations in times to failure across the higher temperature 
range. 

The effects of oxidation were also determined from analysis of polished cross-sections 
and fracture surfaces using optical and scanning electron microscopes. Oxidation was observed to 
occur preferentially along microcracks. Polished cross-sections of samples tested in the lower 
temperature range showed oxidation was relatively minimal and occurred in a rough, uneven 
pattern near the edges, but oxidation was also observed deep into the interior of the section, 
suggesting slow oxidation kinetics. Sections from samples tested in the higher temperature range 
showed oxidation occurring primarily along the edges in a uniform manner so that a reaction front 
could be seen. This gave a shrinking core type of effect as outer fibers were consumed as the 
oxidation process moved inward. This pattern along with the observation that fibers deep in the 
interior of the sample remained free of attack suggests the oxidation kinetics are much faster than 
in the lower temperature range. This change in the rate controlling mechanism from the lower 
temperature regime to the higher temperature regime has been attributed to a shift from the kinetics 
being reaction rate limited to being diffusion limited. These experimental results were used to 
provide a basis for developing a model so that a better understanding of diffusion controlled and 
reaction controlled kinetics could be gained. 

A finite difference model was developed which simulates the diffusion of oxygen into a 
matrix crack bridged by an array of carbon fibers. The model provides a visual, qualitative and 
quantitative basis for understanding the oxidation kinetics. The influence of important variables: 
temperature (T), diffusion coefficient (D), and reaction rate constant (K), are studied. 

By obtaining a better understanding of the oxidation kinetics of carbon fibers in a ceramic 
matrix, the best approaches for protecting against oxidation can be determined depending on the 
application conditions: environment, temperature, stresses and required component lives. One 
possible approach for protecting the fibers might involve applying a less reactive interphase to the 
fiber. Another may be to include an oxidation inhibitor in the matrix that will react with oxygen to 
form oxides or glasses to seal cracks or to consume oxygen to provide more time before the carbon 
is attacked. Still another approach could be to identify external seal coatings that can better protect 
the interior of the material. One approach may be preferred over another due to exposure 
conditions, cost and/or mission requirements. 


THEORY 


The related theory used in developing the finite difference, cracked matrix model will be 
briefly discussed. In a ceramic matrix composite, oxygen diffuses down a pore or crack and reacts 
with carbon to form an oxide (either carbon monoxide or carbon dioxide). The oxide then diffuses 
out of the pore to the atmosphere. Non-Reactive Matrix 
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Figure 1. 
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Oxygen diffusing into a pore and reacting with carbon. 
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The oxidation kinetics for the recession of carbon can be described by the linear- 
parabolic rate law: 

+ x = t, (1) 

kp k] 

where x (m) is the recession distance of the carbon, k p is the parabolic rate constant (m 2 /sec), k, is 
the linear rate constant (m/sec) and t is time (sec). 

The parabolic and linear rate constants can be related to more familiar mass transport 
constants, the diffusion coefficient and the reaction rate constant respectively. Only the relevant 
equations will be discussed. More thorough discussions and derivations can be found in papers by 
Eckel et al. [7] and Glime et al. [8], 

When the linear term in equation ( 1 ) dominates and carbon dioxide is the evolving oxide, 
the equation for the recession distance is expressed as [7]: 

x = t 1/2 k p 1/2 = t 1/2 (2D C ox / p c ) such that k p = (2D C ox / p c ) (2) 

where D is the diffusion coefficient of oxygen in the pore (m 2 / sec), C ox is the oxygen 
concentration of the atmosphere (mol / m 3 ) and p c is the molar density of carbon (mol / m 3 ). 

In the case where carbon monoxide is the reaction product, Equation (2) for the recession 
distance will have a different form. However, the particular oxide that forms and the exact value 
for the diffusion coefficient is not of critical importance hi this discussion and in the model. The 
main concern at this point in the development of the model is to identify trends. The discussion 
does not consider Knudsen effects in which molecular- wall collisions are more frequent than inter- 
molecular collisions. Knudsen effects are not a factor until the pore becomes as small as 0. 1 pm 
[7], In the material studied, typical cracks in the matrix and seal coating are on the order of 1 pm. 

The second limiting case occurs when k p is much larger than kj. The recession distance 
can be expressed as [8]: 

x = t ki = t K C ox / p c such thatki = K C ox / p c (3). 

The above equation takes an alternate form when the oxygen concentration is substituted 
by the oxygen partial pressure, %, times the total concentration of gas molecules in the atmosphere, 
C T (mol/m 3 ). The total concentration can be substituted using the ideal gas law so that Equation 
(3) becomes [7]: 

k] = (1/ Pc )(X C T ) (K) = (1/ p c ) ( X P/RT) K (4) 

where P is pressure (Pa), R is the gas constant (J/mol K), and T is the absolute temperature (K). 
Substituting values for k p and k] , Equation ( 1 ) becomes 

t = 2^ + x = x 2 + x (5). 

k p k, (2D C ox / p c ) (1/pcX XP/RT)K 

The diffusion coefficient, D, can be obtained from Chapman-Enskog kinetic theory [9] and the 
reaction rate constant, K, can be expressed in an Arrhenius form: 
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(6) & (7) 


= (5.9543 x 10 24 ) T 3/2 (l / + 1 / MM lg , K = k () exp(-Q/RT) 

P cr ab £2ab 

where Dab (m 2 / sec) is the diffusion coefficient of gas A (oxygen) diffusing in the stationary or 
stagnant atmosphere of gas B (oxide), T is temperature (K), M A and M B are the molecular weights 
of gases A and B (kg/mol), and P is pressure (Pa). The symbols ct A b (m) and Dab (dimensionless) 
are Leonard Jones potentials for the collision diameter and collision integral respectively . Q is the 
activation energy (J/mol). 

From the above equations for the diffusion coefficient (6) and the reaction rate constant 
(7), it is seen that one of the two arguments in Equation (5) will dominate depending on 
temperature. The diffusion coefficient depends on T 3/2 , with the effect being even less when the 
decreasing molar density of the gas is considered with increasing temperature. The reaction rate 
constant varies exponentially with temperature when large activation energies are considered. It 
can be seen that Equation (5) describes a two step process occurring in series. Each step has the 
potential to limit the overall process. When the reaction rate constant is small and the diffusion 
coefficient is large, the process will be limited by the reaction rate, i.e. reaction controlled. 

Oxygen will diffuse into the pore faster than it can react so that the pore becomes saturated in 
oxygen. However, if the diffusion coefficient is small compared to the reaction rate constant, 
carbon/oxygen reactions occur as soon as oxygen is supplied. In this case the process is controlled 
by the slower diffusion step, i.e. diffusion controlled. 


MODELING 

In order to gain a better understanding of the diffusion and reaction kinetics observed in 
the experiments, a finite difference model was developed. The physical setup of the model is an 
array of carbon fibers bridging a matrix crack, as illustrated in Figure 2. Figure 3 shows an 
example of an experimental observation of oxidation occurring along fibers bridging a crack, 
which is the physical situation on which the model is based. 




Figure 3. Oxidation of carbon fibers that bridge 
matrix cracks. 
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A 2-Dimensional grid pattern was laid out as a 250x250 mesh of grids, which represents 
an open plane within the bridged space parallel to the two crack surfaces of the matrix. Grids 
within the mesh are designated as either open space containing the diffusing gases, or as carbon 
containing grids. Carbon fibers are represented as a 12x12 array of fibers which have diameter of 
10 grids. Fibers are spaced a distance of 10 grids from one another and from the outer edge. One 
quarter of the grid pattern is shown in Figure 4. 



x 


Figure 4. One quarter plot of the matrix pattern (right) orientated as shown on the left. 

In the model, the outer boundary has a initial oxygen concentration arbitrarily set at 1 .0 
mol/cm 3 , which is assumed to stay constant due to the atmospheric conditions. Oxygen diffuses 
into the matrix. Concentrations are calculated throughout the matrix based on oxidation kinetics in 
which both reactivity and diffusion are both dependent on temperature and one may be limiting the 
overall process. Plots of oxygen concentration as a function of grid position showing gradients 
high or low in oxygen concentration will give clues to the limiting step in the two step process. 

When mass transport occurs by convection, the flux is given by, 

Flux = J = K (C A>1 -C m ) 

where C A1 (which has a value of 0 mol/cm 3 ) is the oxygen concentration at the reaction front or of 
the carbon edge (carbon containing grid) and C „ is the oxygen concentration in the bulk gas 
(oxygen containing grid). 

When mass transport occurs by diffusion, the flux is given by, 


Flux = J = D (C 2 - Q ) / Ax 

where (C 2 - C, ) is the difference in oxygen concentration within the stagnant gas across some 
distance, Ax. An important unitless parameter, the Sherwood number, Sh, is obtained from the 
ratio of the fluxes per the driving force due to concentration differences, 

Sh = F(convection) /AC / F(diffusion) /AC = K Ax / D 

The Sherwood number will be used extensively throughout the model. For low Sh values, low 
temperatures, the kinetics are reaction controlled, and for high Sh values, high temperatures, the 
kinetics are diffusion controlled. 
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Detail of the grid pattern showing half of a carbon fiber and open pore space is illustrated 
in Figure 5. The oxygen concentration of a grid is determined by sampling from its four neighbors. 
The fluxes between the central grid and its 4 neighbors are added. The net flux into and out of the 
central grid is zero. The flux equation is then 
solved for the concentration of the central 
grid. When all four neighbors are within the 
stagnant gas as for grid 1 , mass transport 
occurs by diffusion and the equation for the 
oxygen concentration of this grid is 


C(i,j) = sum of oxygen concentrations (i,j-l) 

4 

When one of the neighboring grids contains (i+Fj) 

carbon, as for grid 2 in Figure 5, mass 
transport occurs through convection between 
the central grid and the carbon containing 
grid. For the other three grids mass transport 
occurs by diffusion. The equation for the 
oxygen concentration of the central grid becomes Figure 5 . Grid pattern for model. 

c(i,j) = iC(i.i-n+C(i-i.i)+qi.i+i)i + sh*cq+i.i) 

(3+Sh) 

When two neighbors contain carbon and two contain the stagnant gas, as in grid 3, the equation 
becomes 

ca,j) = ic(i.i-n+c(i-i.m + shiQi,i+n+c(i+i.i)i 

(2+2Sh) 

The program is run for several different Sh values: 0.0001, 0.001, 0.01, 0.1, 1, 10. 
Oxygen concentrations are determined throughout the bridged section based on the boundary 
condition and the equations for oxygen concentration. The program runs until a quasi-steady state 
is reached. Note the program does not remove carbon by reaction at this point; future development 
of the program will remove carbon and continuously recalculate oxygen concentrations so that 
oxidation patterns relating to diffusion controlled and reaction controlled kinetics are obtained. 

The program may be thought of at this point as having a continuous supply of oxygen and carbon, 
where oxygen concentrations throughout the section are calculated until steady state is reached. It 
should be pointed out that since oxidation really occurs in 3-dimensions, the 2- dimensional model 
does not consider or allow the analysis of certain factors. These include 3rd dimensional effects 
such as fiber recession and the supply of oxygen through open channels where the carbon fiber has 
been consumed. 

The obtained data files for a given Sherwood number were used to prepare plots of the 
oxygen concentrations across the interior of the fiber bridged matrix. Gradients in oxygen 
concentration across the bridged section show regions where oxygen supply is low or high and 
give indications about the reactivity. Figure 6 shows the oxygen concentrations for one-quarter of 
the whole section for a relatively high Sherwood number, Sh=0. 1 . The 3-D effect of the plot is due 
to plotting oxygen concentrations across the 2-D surface. It is seen that oxygen concentrations 
drop sharply upon diffusion into the matrix. This suggests that the carbon is very reactive and 
reacts with oxygen as soon as it is supplied. The overall oxidation process is limited by the slow 
supply of oxygen to the reaction process at the carbon/oxygen interface, i.e. diffusion controlled 
kinetics. The case for a relatively low Sherwood number, Sh=0.001, is shown hi Figure 7. The 
plot shows that the reactions are slow so that oxygen is able to diffuse into the interior hi high 
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concentrations and saturate the interior. In this case, the two step process of oxygen diffusion and 
carbon/oxygen reactions is limited by the slower reaction process, reaction controlled kinetics. 
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Figure 6 (left) Sh=0. 1. Figure 7 (right) Sh=0.001. One -quarter plots of a 250x250 matrix with a 12x12 fiber 
array. Interated to equilibrium. 

In the two previous figures, the curves at the edges of the plots represent the oxygen 
concentration through one-half of the thickness of the section. Full plots of the oxygen 
concentrations across the midline are shown in Figure 8 for several Sherwood numbers. High Sh 
values illustrate how the process is diffusion controlled. Oxygen is quickly consumed before it can 
enter into the interior of the matrix. Low Sh values illustrate how the kinetics are reaction 
controlled. Carbon/oxygen reactions are slow so that there is a high concentration of oxygen 
throughout the section. 


Sh = Ax K 
D 

Reaction Controlled 


Sh = 1 

and 10 | 

"(lower) 

Diffusion Controlled 


0 50 100 150 200 250 

Position Across Matrix (0-250) 

Figure 8. Oxygen concentrations across the center of the surface for several Sherwood (Sh) values. 
SUMMARY 

Experimental studies on C/SiC composite material exposed to stressed oxidation 
conditions suggested two types of oxidation kinetics. Lives of the samples were very temperature 
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dependent in the range of 350°C-750°C, while there was less of a temperature dependence hi the 
range of 1000°C-1500°C. Microstructural analysis revealed a different type of oxidation pattern 
for each temperature range. The related theories for carbon recession, reaction kinetics, diffusion 
kinetics, flux due to mass transport, and Sherwood numbers were reviewed. Temperature effects 
on the reaction rate constant and the diffusion coefficient were also discussed. Based on 
experimental observation, and theory, a finite difference model was presented. 

CONCLUSIONS 

The oxidation of carbon fibers in a ceramic matrix is a two step process with a varying 
temperature dependence. The overall process will be limited by the slower step. The kinetics can 
be analyzed by studying the diffusion of oxygen and determining its concentration throughout the 
section. The finite difference model is in good agreement with experimental results in suggesting 
two types of kinetics. The model supports the conclusion that the kinetics are reaction controlled 
at low temperatures, and diffusion controlled at high temperatures. 

The next steps in the model development include the removal of carbon as it reacts with 
oxygen. In the case of diffusion controlled kinetics, i.e. high Sherwood number, oxygen will 
oxidize away the outer fibers first, then continue to move inward so that a shrinking core effect is 
seen. For reaction controlled kinetics, i.e. low Sherwood number, there will be a generalized 
oxidation of fibers throughout the section. Another developmental step for the model will include 
putting an element around the fibers to represent the interphase. Different reactivities of the fiber 
and interphase can then be considered. 
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